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ABSTRACT 

Context. The growth processes from protoplanetary dust to planetesimals are not fully understood. Laboratory experi- 
ments and theoretical models have shown that collisions among the dust aggregates can lead to sticking, bouncing, and 
fragmentation. However, no systematic study on the coUisional outcome of protoplanetary dust has been performed so 
far so that a physical model of the dust evolution in protoplanetary disks is still missing. 

Aims. We intend to map the parameter space for the collisional interaction of arbitrarily porous dust aggregates. This 
parameter space encompasses the dust-aggrcgato masses, their porosities and the collision velocity. With such a com- 
plete mapping of the collisional outcomes of protoplanetary dust aggregates, it will be possible to follow the collisional 
evolution of dust in a protoplanetary disk environment. 

Methods. We use literature data, perform own laboratory experiments, and apply simple physical models to get a com- 
plete picture of the collisional interaction of protoplanetary dust aggregates. 

Results. In our study, we found four different kinds of sticking, two kinds of bouncing, and three kinds of fragmentation 
as possible outcomes in collisions among protoplanetary dust aggregates. Our best collision model distinguishes between 
porous and compact dust. We also differentiate between collisions among similar-sized and different-sized bodies. All in 
all, eight combinations of porosity and mass ratio can be discerned. For each of these cases, we present a complete colli- 
sion model for dust-aggregate masses between 10^^^ and 10^ g and collision velocities in the range cm s 
for arbitrary porosities. This model comprises the collisional outcome, the mass(es) of the resulting aggregatc(s) and 
their porosities. 

Conclusions. We present the first complete collision model for protoplanetary dust. This collision model can be used for 
the determination of the dust-growth rate in protoplanetary disks. 
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1. Introduction 

The first stage of protoplanetary growth has still not been 
fully understood. Although our empirical knowledge on the 
collisional properties of dust aggregates has considerably 

widened over the past years (Blum & Wurm 2008), there 
is no self- consistent model for the growth of macroscopic 
(lust aggregates in protoplanetary disks (PPDs). A reason 
for such a lack of understanding is the complexity in the col- 
lisional physics of dust aggregates. Earlier assumptions of 
perfect sticking have been experimentally proven false for 
most of the size and velocity ranges under consideration. 
Recent work also showed that fragmentation and porosity 
play important roles in mutual collisions between proto- 
planetary dust aggregates. In their review paper, Blum & 



* This paper is dedicated to the memory of our dear friend and 
coUeague Prithjof Brauer (14th March 1980 - 19th September 
2009) who developed powerful models of dust coagulation and 
fragmentation, and thereby studied the formation of planetesi- 
mals beyond the meter size barrier in his PhD thesis. Rest in 
peace, Frithjof. 



Wurm (2008) show the complex diversity that is inherent 
to the collisional interaction of dust aggregates consisting 
of micrometer-sized (silicate) particles. This complexity is 

the reason why the outcome of the collisional evolution in 
PPDs is still unclear and why no 'grand' tlieory on the for- 
mation of planetesimals, based on firm physical principles, 
has so far been developed. 

The theoretical understanding of the physics of dust ag- 
gregate collisions has seen major progress in recent decades. 
The behavior of aggregate collisions at low collisional en- 
ergies - where the aggregates show a fractal nature - is 
theoretically described by molecular dynamics simulations 
of Dominik & Tielens (1997). The predictions of this model 
- concerning aggregate sticking, compaction, and catas- 
trophic disruption - could be quantitatively confirmed by 
laboratory collision experiments of Blum & Wurm (2000). 
Also, the collision behavior of macroscopic dust aggre- 
gates was successfully modeled by a smooth particle hy- 
drodynamics method, calibrated by laboratory experiments 
(Giittler et al. 2009; Geretshauser et al. 2009). These sim- 
ulations were able to reproduce bouncing collisions, which 
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were observed in many laboratory experiments (Blum & 
Wurm 2008). 

However, as laboratory experiments have shown, colli- 
sions between dust aggregates at intermediate energies and 
sizes are characterized by a plethora of outcomes: ranging 
from (partial) sticking, bouncing, mass transfer, to catas- 
trophic fragmentation (see Blum & Wurm 2008). From this 
complexity, it is clear that the construction of a simple the- 
oretical model that agrees with all these observational con- 
straints is very challenging. However, in order to under- 
stand the formation of planetesimals, it is imperative to 
describe the entire phase-space of interest, i.e., to consider 
a wide range of aggregate masses, aggregate porosities, and 
collision velocities. Likewise, the collisional outcome is a key 
ingredient of any model that computes the time evolution 
of the dust size distribution. These collisional outcomes arc 
mainly determined by the collision velocities of the dust 
aggregates, and these depend on the disk model, i.e. the 
gas and material density in the disk and the degree of tur- 
bulence. Thus, the choice of the disk model (including its 
evolution) is another major ingredient for dust evolution 
models. 

These concerns lay behind the approach we adopt in this 
and subsequent papers. That is, instead of first 'funneling' 
the experimental results through a (perhaps ill-conceived) 
theoretical collision model and then to calculate the colli- 
sional evolution, we will directly use the experimental re- 
sults as input for the collisional evolution model. The draw- 
back of such an approach is of course that experiments on 
dust aggregate collisions do not cover the whole param- 
eter space and therefore need to be extrapolated by or- 
ders of magnitude, based on simple physical models which 
accuracy might be challenged. However, we feel that this 
drawback is more than justified by the prospects that our 
new approach will provide: through a direct mapping of the 
laboratory experiments, collisional evolution models can in- 
crease enormously in their level of realism. 

In Paper I, we will classify all existing dust-aggregate 
collision experiments for silicate dust, including three addi- 
tional original experiments not published before, according 
to the above parameters (Sect. 2). We will show that we 
have to distinguish between nine different kinds of colli- 
sional outcomes, which we physically describe in Sect. 3. 
For the later use in a growth model, we will sort these 
into a mass- velocity parameter space and find that we have 
to distinguish between eight regimes of porous and com- 
pact dust-aggregate projectiles and targets. We will present 
our collision model in Sect. 4 and the consequences for the 
porosities of the dust aggregates in Sect. 5. In Sect. 6, we 
conclude our work and give a critical review on our model 
and the involved necessary simplifications and extrapola- 
tions. 

In Paper II (Zsom et al. 2009), we will then, based 
upon the results presented here, follow the dust evolution 
using a recently invented Monte-Carlo approach (Zsom & 
DuUemond 2008) for three different disk models. This is the 
first fully self-consistent growth simulation for PPDs. The 
results presented in Paper II represent the state-of-the-art 
modeling and will give us important insight into questions, 
such as if the meter-size barrier can be overcome and what 
the maximum dust- aggregate size in PPDs is. i.e. whether 
pebbles, boulders, or planetesimals can be formed. 



2. Collision Experiments with Relevance to 
Planetesimal Formation 

In the past years, numerous laboratory and space exper- 
iments on the colUsional evolution of protoplanetary dust 
have been performed (Blum & Wurm 2008). Here, we con- 
centrate on the dust evolution around a distance of 1 AU 
from the solar-type central star where the ambient temper- 
ature is such that the dominating material class are the 
silicates. This choice of 1 AU reflects the kind of laboratory 
experiments that are included in this paper, which were all 
performed with Si02 grains or other refractory materials. 
The solid material in the outer solar nebula is dominated 
by ices, which possibly have very different material proper- 
ties than silicates, but only a small fraction of laboratory 
experiments have dealt with these colder (ices, organic ma- 
terials) or also warmer regions (oxides). In Sect. 6.2, we 
will discuss the effect that another choice of material might 
potentially have, but as we arc far away from even basi- 
cally comprehending the collisional behavior of aggregates 
consisting of these materials, we concentrate in this study 
on the conditions relevant in the inner solar nebula around 
1 AU. 

Table 1 lists all relevant experiments that address col- 
lisions between dust aggregates of different masses, mass 
ratios, and porosities, consisting of micrometer-sized sili- 
cate dust grains, in the relevant range of collision velocities. 
Experiments 1-16 are taken from the literature (cited in 
Table 1), whereas experiments 17 - 19 are new ones not 
published before. In the following two subsections we will 
first review the previously published experiments (Sect. 2.1) 
and then introduce the experimental setup and results of 
new experiments that were performed to fill some regions 
of interest (Sect. 2.2). All these collisions show a diversity 
of different outcomes for which we classify nine different 
collisional outcomes as displayed in Fig. 1. Details on these 
collisional outcomes are presented in Sect. 3. 

2.1. A Short Review on Collision Experiments 

We briefiy review published results of dust-collision experi- 
ments here since these determine the collisional mapping in 
Sect. 3 and 4. The interested reader is referred to the review 
by Blum fc Wurm (2008) for more information. All experi- 
ments are compiled and referenced in Table 1 where we also 
list the collision velocities and projectile masses as these 
will be used in Sect. 4. Most of the experiments in Table 
1 (exception: Exp 10) were performed under low gas pres- 
sure conditions to match the situation in PPDs and most of 
the experiments were carried out in the absence of gravity 
(i.e. free falling aggregates or micro-gravity facilities), see 
column 4 of Table 1. For the majority of the experiments, 
spherical monodisperse Si02 monomers with diameters be- 
tween 1.0 /im and 1.9 /im were used; some experiments used 
irregular Si02 grains with a wider size distribution centered 
around ~ 1.0 /zm, and Exp 5 used irregular ZrSi04 with 
monomer diameters in the range 0.2 . . . 1.0 /xm. 

Exp 1 A well-known growth mechanism for small 
dust aggregates is the hit-and-stick growth, in which the ag- 
gregates c;olli(le with such a small kinetic energy that they 
stick at each other upon first contact without any restruc- 
turing. The first experiments to unambiguously show that 
the hit-and-stick process is relevant to protoplanetary dust 
aggregation were those by Wurm & Blum (1998), Blum 
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et al. (1998, 2000, 2002) and Krause & Blum (2004). These 
proved that, as long as the collision velocities for small 
dust aggregates stay well below 100 cm s~^, sticking col- 
lisions lead to the formation of fractal aggregates. This is 
in agreement with the molecular-dynamics simulations by 
Dominik & Tielens (1997) and Wada et al. (2007, 2008, 
2009). The various experimental approaches for Exp 1 ~ 3 
used all known sources for relative grain velocities in PPDs, 
i.e. Brownian motion (Exp 3), relative sedimentation (Exp 
1), and gas turbulence (Exp 2). In these papers it was also 
shown that the hit-and-stick growth regime leads to a quasi- 
monodispcrsc evolution of the mean aggregate masses, de- 
pleting small grains efficiently and rapidly. For collisions be- 
tween these fractal aggregates and a solid or dusty target. 
Blum & Wurm (2000, Exp 4) found growth at even higher 
velocities, in which the aggregates were restructured. This 
is also in agreement with molecular-dynamics simulations 
(Dominik & Tielens 1997), and so this first stage of pro- 
toplanetary dust growth has so far been the only one that 
could be fully modeled. 

Exp 5: Blum & Miinch (1993) performed collision exper- 
iments between free falling ZrSi04 aggregates of intermedi- 
ate porosity {<j) = 0.35, where 4> is the volume fraction of the 
solid material) at velocities in the range 15 - 390 c;ni s"^. 
They never found sticking, but, depending on the collision 
velocity, the aggregates bounced {v < 100 cm s~^) or frag- 
mented into a power-law size distribution (w > 100 cm s"-'^). 
The aggregate masses were varied over a wide range (10~® 
to 7 X 10~^ g) and also the mass ratio of the two collision 
partners ranged from 1:1 to 1:66. The major difference to 
experiments 1-4 which inhibited sticking in these colli- 
sions were the aggregate masses and their non-fractal but 
still very porous nature. 



Exp 6 - 8: A new way of producing highly porous, 
macroscopic dust aggregates (</> = 0.15 for 1.5 /xm diam- 
eter Si02 monospheres) as described by Blum & Schrapler 
(2004) allowed new experiments, using the 2.5 cm diameter 
aggregates as targets and fragments of these as projectiles 
(Langkowski et al. 2008, Exp 6). In their collision experi- 
ments in the Bremen drop tower, Langkowski et al. (2008) 
found that the projectile may either bounce off from the 
target at intermediate velocities (50 - 250 cm s~^) and ag- 
gregate sizes (0.5 - 2 mm), or stick to the target for higher 
or lower sizes and velocities, respectively. This bouncing 
went with a previous slight intrusion and a mass transfer 
from the target to the projectile. In the case of small and 
slow projectiles, the projectile stuck to the target, while 
large and fast projectiles penetrated into the target and 
were geometrically embedded. They also found that the 
surface roughness plays an important role for the sticking 
efficiency. If a projectile hits into a surface depression it 
sticks while it bounces off when hitting onto a hill with a 
small radius of curvature comparable to that of the projec- 
tile. A similar behavior for the sticking by deep penetration 
was also found by Blum & Wurm (2008, Exp 7) when the 
projectile aggregate is solid - a mm-sized glass bead in their 
case. Continuous experiments on the penetration of a solid 
projectile (1 to 3 mm diameter) into the highly porous tar- 
get {(j) = 0.15, Blum & Schrapler 2004) were performed by 
Giittler et al. (2009, Exp 8) who studied this setup for the 
calibration of a smoothed particle hydrodynamics (SPH) 
collision model. We will use their measurement of the pen- 
etration depth of the projectile. 

Exp 9 - 10: As a. follow-up experiments of the study of 
Blum fc Miinch (1993), D. Heifielmann, H.J. Fraser and J. 

Blum (in prep., Exp 9) used 5 mm cubes of these highly 
porous (0 = 0.15) dust aggregates and colUded them with 
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Fig. 1. We classify the variety of laboratory experiments into nine kinds of coUisional outcomes, involving sticking (S), bouncing 
(B) and fragmenting (F) collisions. All these coUisional outcomes have been observed in laboratory experiments and detailed 
quantities on the outcomes can be given in Sect. 3. These experiments also showed that we have to distinguish between collisions 
of similar-sized aggregates (left columns) and different-sized aggregates (right columns) and some kinds only occur for one of these 
cases (e.g. S3). 



each other {v = 40 cm s~^) or with a compact, cp = 0.24, 

dust target (?; = 20 cm s~^). In both cases they also found 
bouncing of the aggregates and were able to confirm the 
low coefficient of restitution (tiaftcr/^'bciorc) of e = 0.2 for 
central collisions. In their experiments they could not see 
any deformation of the aggregates, due to the limited reso- 
lution of their camera, which could have explained the dis- 
sipation of energy. This was followed up by Weidling et al. 
(2009, Exp 10) who studied the compaction of the same 
aggregates, which repeatedly collided with a solid target. 
They found that the aggregates decreased in size (without 
losing significant amounts of mass) which is a direct mea- 
surement of their porosity. After only 1 000 collisions the 
aggregates were compacted by a factor of two in volume 
filling factor and the maximum filling factor for the veloc- 
ity used in their experiments (1 30 cm s^^) was found 
to be = 0.36. In four out of 18 experiments, the aggre- 
gate broke into few pieces and they derived a fragmentation 
probability of Pfrag = 10~* for the aggregate to break in a 
collision. 

Exp 11: Also using fragments of the high porosity 
{4> = 0.15) dust aggregates of Blum & Schrapler (2004) 
as well as intermediate porosity (0 = 0.35) aggregates, 
Lammel (2008, Exp 11) followed up the fragmentation ex- 



periments of Blum & Miinch (1993). For velocities from 320 

to 570 cm he found fragmentation and measured the 
size of the largest fragment as a measure for the fragmen- 
tation strength. 

Exp 12 - 14-: Exposing the same highly porous (0 = 
0.15) dust aggregate to a stream of single monomers with 
a velocity from 1500 to 6 000 cm s~^, R. Schrapler and J. 
Blum (in prep., Exp 12) found a significant erosion of the 
aggregate. One monomer impact can easily kick out tens 
of monomers for the higher velocities examined. From an 
analytic model, they estimated the minimum velocity for 
this process to be approx. 350 cm s^^. On a larger scale. 
Wurm et al. (2005a, Exp 13) and Paraskov et al. (2007, 
Exp 14) impacted dust projectiles with masses of 0.2 to 
0.3 g and solid spheres into loosely packed dust targets. In 
the drop-tower experiments of Paraskov et al. (2007) they 
were able to measure the mass loss of the target which was - 
velocity dependent - up to 35 projectile masses. The lowest 
velocity in these experiment was 350 cm s~^. 

Exp 15 - 16: In a collision between a projectile of inter- 
mediate porosity and a compressed dust target at a velocity 
above 600 cm s~^, Wurm et al. (2005b, Exp 15) found frag- 
mentation of the projectile but also an accretion of mass 
onto the target. Depending on the collision velocity, this 



Guttler at al.: The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals? 



5 



accretion was up to 0.6 projectile masses in a single colli- 
sion. Teiser & Wurm (2009a, Exp 16) studied this partial 
sticking in many collisions, where solid targets of variable 
sizes were exposed to 100 to 500 /im diameter dust aggre- 
gates with a mean velocity of 770 cm s^^. Although they 
cannot give an accretion efficiency in a single collision, they 
found a large amount of mass accretion onto the targets, 
which is a combination of the pure partial sticking and the 
effects of the Earth's gravity. Teiser & Wurm (2009a) argue 
that this acceleration is equivalent to the acceleration that 
micron-sized particles would experience as a result of their 
erosion from a much bigger body that had been (partially) 
decoupled from the gas motion in the solar nebula. 

2.2. New Experiments 

In this section, we will present new experiments which we 
performed to fill some parameter regions where no pub- 
lished data existed so far. All experiments cover collisions 
between porous aggregates with a solid target and were per- 
formed with the same experimental setup, consisting of a 
vacuum chamber (less than 0.1 mbar pressure) with a dust 
accelerator for the porous projectiles and an exchangeable 
target. The accelerator comprises a 50 cm long, 3 cm di- 
ameter plastic rod in a vacuum feed through. The pressure 
difference between the ambient air and the pressure in the 
vacuum chamber drives a constant acceleration, leading to 
a projectile velocity of up to 900 cm s"-'^, when the accel- 
erator is abruptly stopped. The porous projectile flies on 
and collides either with a solid glass plate (Sect. 2.2.1 and 
2.2.2) or with a free falling glass bead which is dropped 
when the projectile is accelerated (Sect. 2.2.3). The col- 
lision is observed with a high-speed camera to determine 
aggregate and fragment sizes and to distinguish between 
the coUisional outcomes (i.e. sticking, bouncing, and frag- 
mentation). The experiments in this section are also listed 
in Table 1 as Exp 17 to 19. 

2.2.1. Fragmentation with Mass Transfer (Exp 17) 

In this experiment, mm-sized aggregates of different vol- 
ume filling factors {(j) = 0.15 and (f> = 0.35) collided with a 
flat and soUd glass target and fragmented as the colUsion 
velocity was above the fragmentation threshold of approx. 
100 cm s~^. The projected projectile size and its veloc- 
ity were measured by a high-speed camera (see Fig. 2). In 
few experiments, the sizes of the produced fragments were 
measured for those fragments that were sharply resolved, 
which yielded a size distribution of a representative number 
of fragments (the number of resolved fragments varied from 
100 to 400). Assuming a spherical shape of the fragments 
and an unchanged porosity from the original projectile, we 
calculated a cumulative mass distribution as shown in Fig. 
3, where the cumulative mass fraction Y^i^Qinii/Mp) is 
plotted over the normalized fragment mass m.^^jm^. Here, 

TTii and Mp = ^Yl!i=i ^^'^ the mass of the i-th smallest 
fragment and the total mass of all visible fragments and 
N is the total number of fragments. We found that the 
cumulative distribution can well be described by a power 
law 

J\{m')m' dm' = (^^^ , (1) 




Fig. 2. Example for a collision of a porous (0 = 0.35) aggregate 
with a solid target at a velocity of 620 cm s~^. The aggregate 
fragments according to a power-law size distribution and some 
mass sticks to the target (bottom frame). 
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Fig. 3. Mass distribution for two experiments at the velocities 
of 120 and 640 cm s~^. For the higher masses, the distribution 
follows a power law while the lower masses are depleted due to 
the finite camera resolution. The slopes are the same for both 
experiments and there is only an offset (pre-factor) between the 
two. The inset describes this pre-factor n (cf. Eq. 1) which is a 
measure for the strength of the fragmentation. The value clearly 
decreases with increasing velocity (Eq. 2). 

where m' and m are the mass of the fragments in units of 
the projectile mass and /it is a parameter to measure the 

strength of fragmentation, being deflned as the mass of the 
largest fragment divided by the mass of the original pro- 
jectile. The deviation between data and power-law for low 
masses (see Fig. 3) is due to the finite resolution of the cam- 
era, which could not detect fragments with sizes <C 50 /zm. 
In the 10 experiments where the mass distribution was de- 
termined, the power-law index k was nearly constant from 
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Fig. 4. Mass gain of a solid target in 133 collisions (S. Kothe, C. 
Guttler & J. Blum, unpublished data). The target was weighed 
after every 19 collisions. After 57 collisions, one projectile mass 
of dust chipped off from the target, which is a clear effect of 
gravity. Thus, we added this mass to the following measurements 
(triangles) and fitted a linear mass gain, which is 0.023 x nip in 
every collision (solid line). 

0.64 to 0.93, showing no dependence in velocity which was 
varied from 120 to 840 cm s^^. However, a clear depen- 
dence on the velocity was found for the parameter /i, which 
decreased with increasing velocity as shown in the inset of 
Fig. 3. This increasing strength of fragmentation can be 
described as 

where the exponent has an error of ±0.2. The curve was 
fitted to agree with the observed fragmentation threshold 
of 100 cm s-^ 

It is important to know that the number density of frag- 
ments of a given mass follows from Eq. 1 as 

n(m') = ^m'^-\ (3) 

and that the power law for this mass distribution can be 
translated into a power-law size distribution n{a) oc with 
A = 3k — 4. This yields A values from —2.1 to —1.2, much 
flatter than the power-law index of —3.5 from the MRN 
distribution (Matliis et al. 1977), which is widely used for 
the description of high-speed fragmentation of solid materi- 
als. Moreover, this power-law index is consistent with mea- 
surements of Blum & Munch (1993) who studied aggregate- 
aggregate collisions between millimeter-sized ZrSi04 aggre- 
gates (see Sect. 2). Their power-law index, eriuivalent to A 
was —1.4, and for different velocities they also found a con- 
stant power-law index and a velocity-dependent pre-factor 
(their Fig. 8a). 

While most of the projectile mass fragmented into a 
power-law distribution, some mass fraction stuck to the tar- 
get (see bottom frame in Fig. 2). Therefore, the mass of the 
target was weighed before the collision and again after 19 
shots on the same spot. The mass of each projectile was 
weighed which yielded a mean value of 3.34 ± 0.84 mg per 




Fig. 5. Examples for the experimental outcomes in the collisions 
of small aggregates with a solid target. The collision can lead 
to sticking, bouncing, or fragmentation (from left to right). The 
time between two exposures is 2 ms. 



projectile. The increasing mass of the target in units of the 
projectile mass is plotted in Fig. 4. After 57 collisions, dust 
chipped off the target which can clearly be accounted to the 
gravitational influence. For the following measurements we 
therefore added one projectile mass to the target because 
we found good agreement with the foregoing values for this 
offset. The measurements were linearly fltted and the slope, 
which determines the mass gain in a single collision, is 2.3 % 
(S. Kothe, C. Giittler & J. Blum, unpublished data). 

2.2.2. Impacts of Small Aggregates (Exp 18) 

Using exactly the same setup as in the previous section, we 
performed collision experiments with very small (20 fxm to 
1.4 mm diameter) but non-fractal projectiles. Those aggre- 
gates were fragments of larger dust samples as described by 
Blum & Schrapler (2004) and had a volume filling factor 
of (6 = 0.15. In this experiment we observed not only frag- 
mentation but also bouncing and sticking of the projectiles 
to the solid glass target. Thus, the analysis with the high- 
speed camera involved the measurement of projectile size, 
collision velocity, and coUisional outcome, where we dis- 
tinguished between (1) perfect sticking, (2) perfect bounc- 
ing without mass transfer, (3) fragmentation with partial 
sticking, and (4) bouncing with partial sticking. The dif- 
ference between the cases (3) and (4) is that in a fragmen- 
tation event at least two rebounding aggregates were pro- 
duced, whereas in the bouncing collision only one aggregate 
bounced off. 

For the broad parameter range in diameter (20 to 
1400 /xm) and velocity (10 to 1000 cm s^^), we per- 
formed 403 individual collisions in which we were able to 
measure size, velocity, and collisional outcome. Examples 
for sticking, bouncing, and fragmentation are shown in 
Fig. 5. The full set of data is plotted in Fig. 6, where 
different symbols were used for different collisional out- 
comes. Clearly, collisions of large aggregates and high ve- 
locities lead to fragmentation, while small aggregates rather 
bounce off the target. For intermediate aggregate mass (i.e. 
nip = 10^^ g), all kinds of collisions can occur. The back- 
ground color shows a sticking probability which was calcu- 
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Fig. 6. Overview on collision experiments between 20 to 
1400 /im diameter aggregates and a solid target, which leads 
to sticking (diamonds), bouncing (triangles), or fragmentation 
(crosses). The intermediate sticking-bouncing collision is indi- 
cated by the squared symbols. The color indicates the sticking 
probability, i.e. the fraction of sticking events in a logarithmic 
bin around every node. The dotted box denotes the approxi- 
mated parameter range and the solid lines denote the threshold 
between sticking, bouncing and fragmentation as also used in 
Fig. 11. 



lated as a boxcar average (logarithmic box) at every node 
where an experiment was performed. Blue color denotes 
a poor sticking probability while a green to yellow color 
shows a sticking probability of approx. 50 %. We draw the 
solid lines in a polygon [(100,70,800,200,200,17) cm s-\ 
(1.6 • 10-4, 5 • 10-^ 1 • 10-^ 8 • 10-1", 1 • 10-^ 1 • 10-8) g] to 
mark the border between sticking and non-sticking as we 
will use it in Sect. 4. For the higher masses, this accounts 
for a bouncing-fragmentation threshold of 100 cm s-^ at 
1.6 • 10-* g (Exp 18) and for the lower masses, we assume 
a constant fragmentation threshold of 200 cm s-^, which is 
in rough agreement with the restructuring-fragmentation 
threshold of Blum & Wurm (2000, Exp 4). For lower veloc- 
ities outside the solid-line polygon, bouncing collisions arc 
expected, whereas for higher velocities outside the poly- 
gon, we expect fragmentation. Thus, an island of enhanced 
sticking probability for 10-^ - 10-'' g aggregates at a broad 
velocity range from 30 to 500 cm s-^ was rather unexpected 
before. The dotted box is just a rough borderline showing 
for which parameters the experiments were performed as it 
will also be used in Sect. 4. 



2.2.3. Collisions Between Similar Sized Solid and Porous 
Aggregates (Exp 19) 

In a collision between a free falling glass bead of 1 mm 
diameter and a porous (0 = 0.15) dust aggregate of 1.5 
to 8.5 mg mass, we observed fragmentation of the porous 
aggregate while some mass was growing on the solid and in- 
destructible glass bead (S. OUiges & J. Blum, unpublished 
data). In this case, the high-speed camera was used with 
a 3D optics that allowed to image the collision from two 
angles, separated by 90°. On the one hand, this made it 
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Fig. 7. The volume gain of a solid particle colliding with a 
porous aggregate depends on the collision velocity. The data 
points are mean values of 11, 8, and 7 individual experiments 
(left to right), thus, the error bars show the la standard devi- 
ation of velocities and volume gain in these. The images with 
a width of 1.9 mm show the original 1 mm glass bead and ex- 
amples for the mass gain in the three corresponding collision 
velocities (S. OUiges & J. Blum, unpublished data). 



possible to exactly measure the impact parameter b, also 
if the offset of the two collision partners is in the line of 
sight of one viewing angle. Moreover, observing the mass 
growth of the solid projectile is not only a projection in 
one direction but can be reconstructed to get a 3D mea- 
surement. So, the relative velocity and aggregate size were 
measured from the images before the collision and the mass 
gain of the solid glass bead was measured after the collision. 
Figure 7 shows a diagram of volume gain in units of projec- 
tile volume (projectile: porous aggregate) over the collision 
velocity. The three data points are averaged over a number 
of experiments at the same velocity. The error bars denote 
the Ifj standard deviation of collision velocities and pro- 
jectile volume, respectively. A clear trend shows that the 
volume gain of the solid particle decreases with velocity 
and we fitted the data points with 



AV^ = 0.59 - 6.3 X 10" 



(4) 



where Vp is the volume of the glass bead. In this experiment 
we were not able to measure the size distribution of the 
fragments because the absolute velocity is determined by 
the projectile velocity (up to 600 cm s^^), and the faster 
fragments were out of the frame before they were clearly 
separated from each other. 

3. Classification of the Laboratory Experiments 

In this section, the experiments outlined above will be cat- 
egorized according to their physical outcomes in the re- 
spective collisions. In Sect. 2, we saw that various kinds of 
sticking, bouncing, and fragmentation can occur. Here, we 
will keep all these experiments in mind and classify them 
according to nine kinds of possible coUisional outcomes that 
were observed in laboratory experiments. These coUisional 
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outcomes arc displayed in Fig. 1. The denomination of the 
classification follows S for sticking, B for bouncing, and F 
for fragmentation. S and F are meant with respect to the 
target, i.e. the more massive of the two collision partners. 
We will discuss each of the pictograms in Fig. 1, describe 
the motivation for the respective collisional outcomes and 
physically quantify the outcome of these collisions. 

(1) Sticking Collisions: A well known growth mechanism 
is due to hit-and-stick (SI) collisions. Hit-and-stick growth 
was observed in the laboratory (Blum & Wurm 2000; Blum 
et al. 2000) and numerically described (Dominik & Tielens 
1997). Experiments show that the mass distribution during 
the initial growth phase is always quasi-monodisperse. The 
evolution of the mean mass within an ensemble of dust ag- 
gregates due to hit-and-stick (Sl)collisions was calculated 
to follow a power law in time, in good agreement with the 
experiments (Wurm & Blum 1998; Krause & Blum 2004). 
Dominik & Tielens (1997) showed theoretically and Blum 
& Wurm (2000) confirmed this experimentally that small 
fractal aggregates stick at first contact if their collision en- 
ergy is smaller than a threshold energy. For higher energies, 
experiments showed that an aggregate is elastically and 
plastically deformed at the contact zone (Blum & Miinch 
1993; Weidling et al. 2009). This increases the number of 
contacts, which then can lead to sticking at higher veloci- 
ties, an effect we call sticking through surface effects (S2). 
Langkowski et al. (2008) found that sticking can occur for 
even larger velocities, if the target aggregate is porous and 
significantly larger than the projectile. In this case, the pro- 
jectile sticks by deep penetration (S3) into the target and 
cannot rebound, simply because of geometrical consider- 
ations. This effect holds also true if the projectile aggre- 
gate is compact, which has been shown by Blum & Wurm 
(2008) and further studied by Guttler et al. (2009). In Sect. 
2.2.1, we saw that the growth of a solid target can occur 
if a porous projectile fragments and partially sticks to the 
target surface (S4). This growth mechanism was already 
described by Wurm et al. (2005b). Teiser & Wurm (2009a) 
found it to be an efficient growth mechanism in multiple 
collisions. 

(2) Bouncing Collisions: If the collision velocity of two 
dust aggregates is too low for fragmentation and too high 
for sticking to occur, the dust aggregates will bounce (Bl). 
D. Heifielmann et al. (in prep.) found highly inelastic 
bouncing between similar-sized porous dust aggregates and 
between a dust aggregate and a dusty but rather compact 
target, where 95 % of the kinetic energy were dissipated. 
Weidling et al. (2009) showed that the energy can effec- 
tively be dissipated by a significant (and for a single col- 
lision undetectable) compaction of the porous aggregates 
after multiple collisions (collisional outcome bouncing with 
compaction (Bl)). Another kind of bouncing occurred in 
the experiments of Langkowski et al. (2008) in which a 
porous projectile collided with a significantly bigger and 
also highly porous target aggregate. If the penetration of 
the aggregate was too shallow for the S3 sticking to occur, 
the projectile bounced off and took away mass from the 
target aggregate. This bouncing with mass transfer (B2) 
was also observed in the case of compact projectiles (Blum 
& Wurm 2008). 

(3) Fragmenting Collisions: Fragmentation (Fl), i.e. the 
breakup of the dust aggregates, occurs in collisions between 
similar-sized dust aggregates at a velocity above the frag- 
mentation threshold. Blum & Miinch (1993) showed that 



both aggregates are then disrupted into a power-law size 
distribution. If a target aggregate is exposed to impacts 
of single monomer grains or very small dust aggregates, R. 
Schrapler & J. Blum (in prep.) found that the target aggre- 
gate is efficiently eroded (F2) if the impact velocities exceed 
1 500 cm s~^. This mass loss of the target was also observed 
in the case of larger projectiles into porous targets (Wurm 
et al. 2005a; Paraskov et al. 2007). Similar to the Fl frag- 
mentation, it may occur that one aggregate is porous while 
the other one is compact. In that case, the porous aggre- 
gate fragments but cannot destroy the compact aggregate. 
The compact aggregate accretes mass from the porous ag- 
gregate (Sect. 2.2.3). We call this fragmentation with mass 
transfer (F3). 

These nine fundamental kinds of collisions are all based 
on firm laboratory results. Future experiments will almost 

certainly modify this picture and potentially add so far un- 
known collisional outcomes to this list. However, at the 
present time this is the complete picture of possible col- 
lisional outcomes. In the following we will quantify the 
thresholds and boundaries between the different collision 
regimes as well as physically characterize the collisional out- 
comes therein. 

SI: Hit-and-Stick Growth 

Hit-and-stick growth occurs when the collisional energy in- 
volved is less than 5 • i^rou (Dominik & Tielens 1997; Blum 
& Wurm 2000), where E'roii is the energy which is dissi- 
pated when one dust grain rolls over another by an angle of 
90°. We can calculate the upper threshold velocity for the 
hit-and-stick mechanism of two dust grains by using the 
definition relation between rolling energy and rolling force, 
i.e. 



-E'roll = 



(5) 



Here, aq is the radius of a dust grain and Fioii is the rolling 
force. Thus, we are inside the hit-and-stick regime if 



(6) 



where is the reduced mass of the aggregates. The hit- 
and-stick velocity range is then given by 



V < 



. Trao-F'roii 



(7) 



S2: Sticking by Surface Effects 

For velocities exceeding the hit-and-stick threshold velocity 
(Eq. 7), we assume sticking because of an increased contact 
area due to surface fiattening and, therefore, an increased 
number of sticking grain-grain contacts. For the calculation 
of the contact area, we take an elastic deformation of the 
aggregate (Hertz 1881) and get a radius for the contact area 
of 



So 



9 2 



G 



(8) 



Here, v is the collision velocity, G is the shear modulus, 
and is the reduced radius. We estimate the shear mod- 
ulus with the shear strength, which follows after Sirono 
(2004) as the geometric mean of the compressive strength 
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and the tensile strength. These parameters were measured 
by Blum & Schrapler (2004) to be 4000 dyn cm'^ (com- 
pressive strength) and 10 000 dyn cm~^ (tensile strength), 
so we take 6 320 dyn cm^^ for the shear modulus, which is 
consistent with estimates of Weidling et al. (2009). 

The energy of a pair of bouncing aggregates after the 
collision is ^ 

Ercst.= e^-m^.v'^ (9) 



with the coefficient of restitution e. The contact energy of 
the flattened surface in contact is 



Er.. 



(10) 



where Eq is the sticking energy of a monomer grain with 
radius oq. We expect sticking for £^cont. > ^^rest., thus. 
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(12) 



This is the sticking threshold velocity for sticking through 
surface effects (S2), which is based on the Hertzian defor- 
mation which is of course a simplified model but has proven 
as a good concept in many attempts to describe slight defor- 
mation of porous dust aggregates (Langkowski et al. 2008; 
Weidling et al. 2009). 

We have to ensure that the centrifugal force of two ro- 
tating aggregates, sticking like above, does not tear them 
apart, which is the case if 



-f'cent > TtTsI, 



(13) 



where T is the tensile strength of the aggregate material. 
The centrifugal force in the worst case of a perfectly grazing 
collision is 



2, ,2 



2a„ 



(14) 



where 2a^ is a conservative estimation for the radial dis- 
tance of the masses with tangential velocity ev. Thus, only 
collisions with velocities 



V < 







3 


■27rTa^" 


(i) 


G 







(15) 



can lead to sticking. For the relevant parameter range (see 
Table 2 below), the threshold velocity in Eq. 15 is always 
signific;antly greater than the sticking velocity in Eq. 12, 
thus, we can take Eq. 12 as the relevant velocity for process 
S2. 

We will use this kind of sticking not only within the 
mass and velocity threshold as defined by Eq. 12 but also 
for collisions where we see sticking which can so far not be 
explained by any model like in experiment 6 or 18. For all 
these cases, we assume the porosity of target and projec- 
tile to be unchanged, disregarding any slight compaction as 
needed for the deformation. One exception is the sticking of 
small, fractal aggregates, which clearly goes together with 
a compaction of the projectile (Dominik & Tielens 1997; 
Blum & Wurm 2000). In these cases we assume a projec- 
tile compaction by a factor of 1.5 in volume filling factor as 
there is no precise measurement on this compaction. 



S3: Sticking by Deep Penetration 

If the target aggregate is much larger than the projectile, 
porous and flat, an impact of a (porous or compact) projec- 
tile results in its penetration into the target. Sticking is in- 
evitable if the penetration of the projectile is deep enough, 
i.e. deeper than one projectile radius. In that case, the pro- 
jectile cannot bounce off the target from geometric con- 
siderations. This was found in experiments of Langkowski 
et al. (2008) in the case of porous projectiles and Blum & 
Wurm (2008) in the case of solid projectiles. The result of 
the collision for penetration depths Dp > Op is that the 
mass of the target is augmented by the mass of the projec- 
tile and the volume of the new aggregate reads 

V = Vt-Tral{Dp-ap) + ^Vp (16) 

= Vt+ - Tra^Dp , (17) 

with Vp and Vt being the volume of the projectile and 

target, respectively. We distinguish between compact and 
porous projectiles and take the experiments of Giittler 
et al. (2009) and Langkowski et al. (2008) for impacts into 
(j) = 0.15 dust aggregates and calculate the sticking thresh- 
old velocities. 

For compact projectiles, we use the linear relation for 
the penetration depth of Giittler et al. (2009) 

nipV 



Dr. 



7- 



where 



and A-n = 



(18) 

are the projec- 



tile mass and cross section, respectively. Although Giittler 
et al. (2009) suggest a power-law relation for the penetra- 

, we choose the 



0.23±0.13,,0.89±0.34 



tion depth, i.e. Dp = ^nv^ 
linear relation in Eq. 18 for simplicity which is also in agree- 
ment with the data within the error bars. For such a lin- 
ear fit, the slope to the data in Giittler et al. (2009) is 



7 



8.3-10 ^ cvc? s g ^. We assume sticking for Dp > ap 



and get sticking due to process S3 in the velocity range 

-1 



V > 



'PQ0p 



(19) 



which only depends on the projectile bulk density po and 

filling factor (j)p and not on projectile radius. 

A porous projectile, colliding with a porous target, 
makes a visible indentation into the target aggregate if the 
kinetic energy is E > E^ain, with a material-dependent min- 
imum energy Emin- The crater volume is then given by 



cm 



(20) 



(see Fig. 15 in Langkowski et al. 2008). Again, from ge- 
ometrical considerations, we assume that sticking occurs 
if the projectile penetrates at least one radius deep, thus, 
Vci. > 0.5Vp, where Vp = ^Tra^ is the volume of the projec- 
tile. Thus, 

3 

E 
Et 



(21) 



1 9 ^ / 1 m 
2"^^^^* (27 



V > 



niE^ 



(22) 



(23) 
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For these velocities, the projectile is inevitably embedded 
into the target aggregate. However, if the impact energy is 
less than Ej^in, the collision will not lead to a penetration 
so that the final condition for sticking of a porous projectile 
according to process S3 is 



V > max 



I2E„ 



m 



2^ 



54: Partial Sticking in Fragmentation Events 

As introduced in Sect. 2.2.1, a fragmenting collision be- 
tween a porous aggregate and a solid target can lead to a 
partial growth of the target. The mass transfer from the 
projectile to the target is typically 2.3 % of the projectile 
mass (Fig. 4) and without better knowledge we assume that 
the transferred mass has a volume filling factor of l,50p. 
The remaining mass of the projectile fragments according 
to the power-law mass distribution given in Eq. 3, with the 
fragmentation strength from Eq. 2. 

For a compact projectile aggregate impacting a com- 
pact target, the threshold velocity for the S4 process is 
V = 100 cm s~^ and thus identical to that of the Fl process. 
The fragmentation strength is given by Eq. 36. 

Bl: Bouncing with Compaction 

In a bouncing collision we find compaction of the two colli- 
sion partners. For similar-sized aggregates, the increase of 
the volume filling factor was formulated by Weidling et al. 
(2009, their Eq. 25) to be 



t)+{(f>,v) >0 



(25) 



with i^{v) = vo ■ (w/20 cm s ^) "^^^ , 4>niax{v) = (j^o + A(j) 
{v/20 cm s-i)^^^ and uq = 850, <^o = 0.15, A(f) = 0.215 for 
u < 50 cm s^^. Here. 4'nia.yL is the saturation of the filling 
factor after many collisions, which follows an exponential 
function with the e- folding width u (Weidling et al. 2009). 
In their experiments, v was the velocity of a porous pro- 
jectile colliding with a solid target (infinite mass). In the 
case of similar-sized colliding aggregates, the velocity would 
be 0.5 • V for each aggregate in a center-of-mass system. 
Therefore, we scale the velocity as 



1 + 



V 

1+^ 



(26) 
(27) 



where Vp {vt) is the center-of-mass velocity of the projectile 

(target). In the case of nip ^ nit we have the situation of 
Weidling et al. (2009) with Vp = v, thus, these velocities 
are chosen to calculate the scaling of u{v) and (j)ma,x{v) for 
projectile and target compaction, respectively. This means 
that a projectile with negligible mass with respect to the 
target cannot compact the target but is only compacted by 
itself, while two aggregates of the same mass are equally 
compacted. 

For 0max(^^), Weidling et al. (2009) gave the above rela- 
tion which is biased by the experimentally used dust sam- 
ples and overestimates the compression for very low veloc- 
ities. Therefore, we propose an alternative scaling relation 
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Fig. 8. The original compressive strength curve measured by 
Giittler et al. (2009) (Eq. 29, solid hno) is biased by the dust 
samples used in the experiments. To describe also the compres- 
sion of dust aggregates with a volume filling factor lower than 
those used by Giittler et al. (2009), we extrapolate the curve 
with a power law (Eq. 30, dashed line) for p < pm. 



for 4>niax{v)- In a collision with velocity v we can calculate 
a dynamic pressure 



(28) 



This pressure is increased by a factor i'{v) as we know from 
the experiments of Weidling et al. (2009) that the contact 
area is very small (factor 1/v of the aggregate surface) 
and that only a very confined volume is compressed. For 
?; = 20 cm s~^ the pressure calculated from Eq. 28 is very 
close to the value given by Weidling et al. (2009). From 
this pressure we calculate the compression from the com- 
pressive strength curve which Giittler et al. (2009) derived 
for collisions: 



f^comp 



(P) 



exp 



(29) 



with (t)i = 0.12, 4)2 = 0.58, A = 0.58, and = 1.3 x 
10** dyn cm~^. This compressive strength curve is also bi- 
ased from the experiments as its lowest value is = 0.12. 
Assuming the saturation part of the compressive strength 
curve to be general, we propose a power law for p < 
with the same slope as in Eq. 29 for 0comp(Pm) which is 
then given by 



+ < 



Pn 



02+^1 2AlnlO 



(30) 



and is able to treat the lowest filling factors and pressures. 
Equations 29 and 30 determine the compression in a con- 
fined volume. Taking into ac;count that after many colli- 
sions only an outer rim of the aggregate is compressed, we 
reduce the compression by a factor = 0.79 to fit the 
(f'niaxiv = 20 cm s^^) = 0.365 experimentally measured by 
Weidling et al. (2009). 
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Conclusively, we calculate the increase of the volume 
filling factor from Eq. 25, where ^max is now provided by 
the dynamical pressure curve as 

<?^max(v) = /c • ?!'comp(Pdyn) , (31) 

where (/)comp is given by Eqs. 29 and 30. For the pressure we 
use Eq. 28 and for the corresponding velocities we use Eqs. 
26 and 27 to calculate projectile and target compression, 
respectively. The maximum compression (pmnKiv), which an 
aggregate can achieve in many collisions at a given velocity, 
is shown in Fig. 8. 

Weidling et al. (2009) found that in this bouncing 
regime, the aggregates can also fragment with a low prob- 
ability. We adopt this fragmentation probability of 

Pfrag = 10-" (32) 

and assume that an aggregate breaks into two similar-sized 
fragments as suggested by their Fig. 5. 

B2: Bouncing with Mass Transfer 

Langkowski et al. (2008) and Blum & Wurm (2008) found, 
that the collision between a projectile (porous or solid) and 

a porous target aggregate can lead to a slight penetration 
of the projectile into the target followed by the bouncing of 
the projectile. This leads to a mass transfer from the target 
to the projectile (see Fig. 7 in Langkowski et al. 2008). 
We assume that the transferred mass is one projectile mass 
(Fig. 8 in Langkowski et al. 2008), thus, 

Amt^p = rrip (33) 

and that the filling factor of the transferred (compacted) 
material is 1.5 times that of the original target material, 
i.e. 

0t^p = 1.5 X (l)t . (34) 

Although the filling factor of the transferred material was 

not measured, we know that the material is significantly 
compacted in the collision (see x-ray micro tomography 
(XRT) analysis of Giittler et al. 2009) so that the above 
assumption seems justified. 



Fl: Fragmentation 

When two similar-sized dust aggregates collide at a velocity 
which is greater than the fragmentation velocity of 

Wfrag = 100 cm S^^, (35) 

they will both be disrupted. Blum & Miinch (1993) found 
fragmentation for mm-sized ZrSi04 dust aggregates with a 
porosity oi (j) = 0.35 at a velocity greater than 100 cm s~^. 
In their experiments, the aggregates fragmented according 
to a power-law size distribution with an exponent of A = 
— 1.4 (see Sect. 2.2.1) which we will use hereafter. The two 
largest fragments together have a mass of ii{v){mp + nit), 
where we can determine /x(w) from the experiments of Blum 
& Miinch (1993, ZrSi04 aggregate collisions with (f) = 0.35) 
and Lammel (2008, Si02 aggregates of different porosities). 
These values are plotted in Fig. 9 and a power-law fit for 
velocities v > 100 cm s~^ 

/ V \ -0-31 

fi(v) = ( A (36) 

' VlOO cm s-i/ ^ ' 

is shown by the solid line, which is again fitted to match the 

fragmentation threshold of 100 cm (cp. Eq. 2). Here, 

the error in the exponent is ±0.02. 




□ Blum and MOnch (1993) 
O Lammel (2008), = 0.15 
A Lammel (2008), js=0.35 



100 



impoct velocity [cm/s] 



1000 



Fig. 9. The impact strength for aggregate-aggregate collision 
also increases for higher velocities (decreasing /n, cp. inset in 
Fig. 3). The fitted power law is given by Eq. 36. 

F2: Erosion 

If a projectile collides with a significantly larger porous tar- 
get aggregate at a sufficiently high impact velocity, the tar- 
get may be eroded. R. Schrapler & J. Blum (in prep.) found 
erosion of porous (</> = 0.15) aggregates which were exposed 
to 1.5 fiia diameter Si02 monomers (mass mo) at veloci- 
ties from 1500 to 6 000 cm s~^. Their numerical model, 
which fits the experimental data very well, predicts an on- 
set of erosion for a velocity of 350 cm s~^. The eroded mass 
grows roughly linear with impact velocity, i.e. 



Am 



mi, 



80 



( ^ ) 

VlOO cm s-V 



(37) 



where Am is the amount of eroded mass and rrip = rriQ is 
the projectile mass. Paraskov et al. (2007) also found mass 
loss of a porous target aggregate for velocities from 350 
to 2150 cm s""'^, although the process involved is widely 
different. They used porous and solid projectiles and their 
results (Fig. 4 in Paraskov et al. 2007) are consistent with 



Am 



15 
20 



100 cm s 



(38) 



which is in agreement with non zero-gravity experiments 
of Wurm et al. (2005a), who estimated a mass loss of 10 
projectile masses for velocities larger than 1650 cm s"-*^. 
Due to the small variation in projectile mass within each 
of the two experiments, we apply a power law in mass and 
merge both experiments to 



Am 



6 

80 



0.092 



, 100 cm s" 

The velocity range for erosion is therefore 
Vfr > 350 cm s~^ 



(39) 



(40) 



and is consistent in both experiments. 

For compact targets, R. Schrapler & J. Blum (in prep.) 
were able to measure the velocity range for erosion at 



Ver > 2 500 cm s 



(41) 
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Due to the nature of the compact target, far less material 
was eroded, i.e. 



Am 



8 



550 V 100 cm s 



nop 

mo 



0.092 



(42) 



Here, we applied the same power-law index as in Eq. 39 due 
to the absence of large-scale experiments in this case. We 
assume a mass distribution of the eroded material according 
to Eq. 2. 

F3: Fragmentation with Mass Transfer 

In Sect. 2.2.3 wc described the vohimc transfer from a 
porous aggregate to a solid sphere (assumed to be repre- 
sentative for a compact aggregate) above the fragmentation 
threshold velocity (see Eq. 4). Without better knowledge, 
we assume that the transferred mass has a volume filling 
factor of 1.5 times that of the porous collision partner (c/ip) 
and cannot exceed the mass of the porous aggregate, thus 

Am = mp(t)1.5(?ip (o.59 - 6.3 x 10"^^^^^^) ' (43) 

where mp(t) is the mass if the porous aggregate which can 
either be projectile or target in our definition, depending 
on its actual mass. For the fragmentation of the porous 
aggregate we assume a power-law distribution following the 
Fl case. If the collision velocity is higher than 940 cm s^^, 
Eq. 43 yields no mass gain for the compact aggregate, thus, 
the mass of the compact aggregate is conserved and only 
the porous aggregate fragments. 



4. Collision Regimes 

In this Section wc intend to build on the physical descrip- 
tions, which wc have derived in the previous Section, and 
(icvclop a complete collision model for determination of 
the coUisional outcome in protoplanetary dust interactions 
(Fig. 1). This means that for each collision that may occur, 
a set of collision parameters will be provided as input for 
a numerical model of the evolution of protoplanetary dust 
(see Paper II). The most crucial parameters that mainly 
determine the fate of the colliding dust aggregates in each 
collision are the respective dust-aggregate masses and their 
relative velocity. 

Moreover, in Sect. 2 and 3, we saw that the porosity 
difference between the two collision partners also has a big 
impact on the collisional outcome. The only difference be- 
tween the outcomes Fl and F3 (and between S3 and S4) 
is that the target aggregate is cither porous or compac;t. 
Thus, we define a critical porosity 4>c to distinguish between 
porous or compact aggregates. This value can only roughly 
be confined between = 0.15 (S3 sticking, clearly an effect 
of porosity, Langkowski et al. 2008) and (f) ~ 0.64 (random 
close packing, clearly compact Torquato et al. 2000) and 
without better knowledge we will choose 0c = 0.4. 

Another important parameter is the mass ratio of the 
collision partners. Again, the sticking by deep penetration 
(S3) occurs for the same set of parameters as the fragmen- 
tation (Fl) and only the critical mass ratio Tm — mt/m-p 
is different. From the work of Blum & Miinch (1993) and 
Langkowski et al. (2008) , we can confine this parameter to 
the range 10 < r,,! < 1 000 and will also treat it in Paper II 
as a free parameter (with fixed values = 10, 100, 1 000). 




^ 

.PC CP ^ 

p ^^^^^ 



Fig. 10. Exporiiricnts suggest that collisions between dust ag- 
gregates with different porosities lead to a different outcome 
than those between aggregates of similar porosity. Thus, our 
model distinguishes between porous and compact aggregates, 
which leads to the displayed four types of collisions ( 'pp ', 'pc ', 
'cp\ 'cc') if the collision partners are not too different in size 
(left). The size ratio of projectile and target aggregate was iden- 
tified as another important parameter and we distinguish be- 
tween similar-sized and different-sized collision partners. Thus, 
in addition to the four collision types on the left, impacts of 
projectiles into much larger targets {'pP', 'pC, 'cP', 'cC; the 
target characterized by a capital letter) can also occur (right). 
The boundary between similar-sized and different-sized aggre- 
gates is given by the critical mass-ratio parameter rm- Collisions 
on the left are restricted to rUp < mt < r^nfUp, collisions on the 
right happen for mt > rmrrip. 



A further parameter, which has an impact on the colli- 
sional outcome, is the impact angle but at this stage, due to 
a lack of information of the actual influence of the impact 
angle on the collisional result, we will treat all collisions as 
central collisions. Experiments by Blum & Miinch (1993), 
Langkowski et al. (2008), or Lammel (2008) indicate rather 
small differences between central and grazing collisions so 
that we feel confident that the error due to this simplifica- 
tion is small. Another parameter, which we also neglect at 
this point due to a lack of experimental data, is the sur- 
face roughness of the aggregates. Langkowski et al. (2008) 
showed its relative importance, but a quantitative treat- 
ment of the surface roughness is currently not possible. 

The binary treatment of the parameters 4>c and r,„ leads 
to Fig. 10 whereafter we have four different porous-compact 
combinations and, if we take into account that the collision 
partners can either be similar-sized or different-sized, we 
have a total of eight collision combinations. We will call 
these 'pp\ 'pP\ 'cc', 'cC, 'cp\ 'cP', 'pc', and 'pC. Here, 
the first small letter denotes the porosity of the projectile 
( 'p ' for porous and 'c ' for compact) and the second letter 
denotes the target porosity which can be either similar-sized 
(small letter) or different-sized (capital letter). Aggregates 
with porosities (j) < <pc are 'porous', those with (p > 4>c are 
'compact'. If the mass of the target aggregate mt < rmrrip, 
we treat the collisions as equal-sized, for mt > T'm'Tip, the 
collisions arc treated as different-sized. 

For each combination depicted in Fig. 10, we have the 
most important parameters (1) projectile mass mp and (2) 
collision velocity v, which then determine the collisional 
outcome. As shown in Fig. 11, we treat each combination 
from Fig. 10 separately and define the collisional outcome as 
a function of projectile mass and collision velocity. For the 
threshold lines and the quantitative collisional outcomes we 
use a set of equations, which were given in Sect. 3. For a 
quantitative analysis and application to PPDs (see Paper 
II), knowledge of the material parameters of the monomer 
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Fig. 11. The resulting collision model as described in this paper. We distinguish between similar-sized (left column) and different- 
sized (right column) collision partners, which are either porous or compact (also see Fig. 10). For each case, the important 
parameters to determine the coUisional outcome are the projectile mass and the collision velocity, collisions within green regions 
can lead to the formation to larger bodies while red regions denote mass loss. Yellow regions are neutral in terms of growth. The 
dashed and dotted boxes show where experiments directly support this model. 



dust grains and dust aggregates is required. In Table 2 we 
list all relevant parameters for 1.5 fim Si02 spheres, for 
which most experimental data are available. However, we 
believe that the data in Table 2 is also relevant for most 
types of micrometer-sized silicate particles. 



The only coUisional outcome, which is the same in all 
regimes, is the hit-and-stick (SI) process, which, due to its 
nature, does not depend on porosity or mass ratio but only 
on mass and collision velocity. Thus, all collision combina- 
tions in Fig. 11 have the same region of sticking behavior for 
a mass-velocity combination smaller than defined by Eq. 7. 
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Table 2. Particle and aggregate material properties used for 
generating Fig. 11. 



symbol value reference 

monomer-grain properties: 

oo 0.75 nm 

mo 3.18 X 10"^^ g 

po 2 g cm"^ 

Eo 2.2 X 10"* erg Blum & Wurm (2000), 

Poppe et al. (2000) 
Froii lO""* dyn Heim et al. (1999) 



aggregate properties: 




e 


0.05 


Blum & Miinch (1993), 
D. Heifielmann et al. (in 
prep.) 


G 


6320 dyn cm"^ 


this work 


T 


10" dyn cm-2 


Blum & Schrapler (2004) 




0.40 


this work 


r-m 


10 - 1 000 


this work 


7 


8.3 X 10"-^ s cm^ g"^ 


Giittler et al. (2009) 


Et 


3.5 X 10'' erg 


Langkowski et al. (2008) 


F, ■ 


3.1 X 10"^ erg 


Langkowski et al. (2008) 


01 


0.12 


Giittler et al. (2009) 


02 


0.58 


Giittler et al. (2009) 


A 


0.58 


Guttler et al. (2009) 


Pm 


1.3 X 10" dyn cm"^ 


Giittler ct al. (2009) 


fc 


0.79 


this work 




850 


Woidling et al. (2009) 


A 


-1.4 


this work 



This parameter region is marked in green because hit-and- 
stick (SI) can in principle lead to the formation of arbitrary 
large aggregates. Marked in yellow are collisional outcomes, 
which do not lead to further growth of the target aggregate, 

but conserve the mass of the target aggregate, which is only 
the case for bouncing with compaction (Bl). For simplic- 
ity, the weak fragmentation probabihty of Pfrag = 10~^ (see 
Sect. 3) has been neglected in the coloring. The red-marked 
regions are parameter sets for which the target aggregate 
loses mass. 

The dashed and dotted boxes in Fig. 11 mark the mass 
and velocity ranges of the experiments from Table 1. In 
Paper II, this plot will help us to see in which parameter 
regions collisions occur and how well they are supported by 
experiments. We will now go through all of the eight plots 
in Fig. 11 and explain the choice for the thresholds between 
the collisional outcomes. 

'pp': In addition to the omnipresent hit-and-stick (SI) 
regime, which is backed by experiments 1 3 in Table 1, col- 
lisions of porous projectiles can also lead to sticking through 
surface effects (S2), whose threshold is determined by Eq. 
12. For higher velocities {v > 100 cm s~^, Eq. 35), frag- 
mentation sets in. Bouncing (Bl) and fragmentation (Fl) 
in this regime are well tested by experiments 5, 9, and 11 
in Table 1. 

'pP': As the projectiles are also porous here, we have 
the same sticking through surface effects (S2) threshold as 
in 'pp '. The same collisional outcome (but with compaction 
of the projectile) was found for collisions of small aggregates 
(Blum & Wurm 2000, experiment 4 in Table 1). Langkowski 
et al. (2008) (experiment 6) found the S2 collisional out- 
come for projectile masses 

mp < 10-4 g. (44) 



thus we have a horizontal upper limit for S2 in the 'pp ' plot 
of Fig. 11. Extrapolation of experiment 6 to large aggregate 
masses 

mp > 0.1 g (45) 

results in bouncing with mass transfer (B2) . A linear inter- 
polation between perfect sticking for mp < 10^* g and per- 
fect bouncing for nip > 0.1 g, justified by the sticking prob- 
abilities shown in Fig. 5 of Langkowski et al. (2008), gives a 
sticking probability for the mass range 10~^ g < rUp < 0.1 g 
(striped region in the 'pP' of Fig. 11) of 

Pstic. = 4log,o(^) (46) 

In Sect. 3 we defined the threshold for sticking by deep 
penetration (S3) by Eqs. 23 and 24, which are prominent in 
the 'pP' plot for high velocities. For even higher velocities, 
we have erosion of the porous aggregate (F2), defined by 
the threshold velocity in Eq. 39 and based on experiments 
12 - 14 in Table 1. 

'cc': Our knowledge about collisions between similar- 
sized, compact dust aggregates is rather limited. Blum & 
Miinch (1993) performed collisions between similar-sized 
aggregates with (p = 0.35. Although this is lower than the 
critical volume filling factor (p^ as defined in Table 2, we 
assume a similar behavior also for aggregates with higher 
porosity. Therefore, without better knowledge, we define a 
fragmentation threshold as in the 'pp ' regime, and take the 
hit-and-stick (SI) threshold for low energies. We omit the 
sticking through surface effects (S2) in this regime because 
of the significantly lower compressibility of the compact ag- 
gregates. 

'cC: Also in this collision regime, the experimental 
background is very limited. For low collision energies we 

assume a hit-and-stick (SI) growth, for higher velocities 
bouncing with compaction (Bl) and, if the fragmentation 
threshold {v > 100 cm s^^, Eq. 35) is exceeded, fragmen- 
tation with mass transfer (S4). Based on experiment 12, 
we have an erosion (F2) limit for velocities higher than 
2 500 cm s-i (Eq. 41). 

'cp'and 'pc': These two cases are almost identical with 
the only difference that the compact aggregate can either 
be the projectile or the target (i.e. slightly lower or higher 
in mass than the target aggregate). However, the mass ra- 
tio of both aggregates is within the critical mass ratio r^. 
Besides the already-discussed cases SI, S2, and Bl, we as- 
sume fragmentation 100 cm s~^ (Eq. 35). Due to the nature 
of the collision between a compact and a porous aggregate, 
only the porous aggregate is able to fragment, whereas the 
compact aggregate stays intact. If the compact aggregate 
is the projectile, the target mass is always reduced, thus 
we have fragmentation with mass transfer (F3) from the 
target to the projectile. If the target is compact, it grows 
by fragmentation with mass transfer (S4), if the velocity 
is less than 940 cm (see Eq. 43). For higher velocities, 
Eq. 43 yields no mass gain and so this region is neutral in 
terms of growth. Collisions at high velocities are confirmed 
by experiment 19 in this regime. 

'cP': While small collision energies lead to hit-and-stick 
(SI) higher energies result in bouncing with mass transfer 
(B2) (Exp. 8, Blum & Wurm 2008). This region is confined 
by the sticking by deep penetration (S3) threshold velocity 
as defined in Eq. 19, based on experiment 7 (Giittler et al. 
2009). At even higher velocities above 350 cm s~^ (Eq. 39), 
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we get erosion of the target aggregate as seen in experiments 
12 - 14. 

'pC: This plot in Fig. 11 looks the most complicated 

but it is supported by a large number of experiments. For 
low collision velocities, we again have hit-and-stick (SI) and 
sticking through surface effects (S2), and a transition to 
bouncing with compaction (Bl) for larger collision energies. 
The existence of the Bl bouncing region has been shown 
in experiments 9 and 10 (D. HciBclmann ct al., in prep.; 
Weidling et al. 2009). For higher velocities and masses 
above 1.6 • 10"** g we assume a fragmentation threshold 
of 100 cm s~^ with mass transfer to the target (S4), as 
seen in experiment 16 (Sect. 2.2.1). For lower masses, the 
odd-shaped box of experiment 18 is a direct input from 
Sect. 2.2.2 (see Fig. 6). In the striped region between Bl 
and S4, we found in experiment 18 a sticking probability of 
^stick = 0.5. For lower masses, experiment 4 showed stick- 
ing through surface effects (S2) with a restructuring (com- 
paction) of the projectile. As in the 'pP' regime, we set 
the threshold for a maximum mass to 8 • 10~^° g, while the 
upper velocity threshold - which must be a transition to a 
fragmentation regime (Blum & Wurm 2000) - is 200 cm 
from experiments 4 and 18. 



5. Porosity Evolution of the Aggregates 

Since the porosity of dust aggregates is a key factor for the 
outcome of dust aggregate collisions (Blum & Wurm 2008), 
it is paramount that coUisional evolution models follow its 
evolution (Ormel et al. 2007, Paper II). Therefore, in this 
section, we want to stress on the evolution of the dust ag- 
gregates' porosities and recapitulate the porosity recipe as 
used in Sect. 3. In this paper we have used the volume fill- 
ing factor as a quantitative value, being defined as the 
volume fraction of material (one minus porosity) . In Paper 
II, we will also use the enlargement parameter as in- 
troduced by Ormel et al. (2007), which is the reciprocal 
quantity ^ = <j)~^ . 

Starting the growth with solid dust grains, we have a 
volume filling factor of 1, which will however rapidly fall due 
to the hit-and-stick (SI) growth, producing highly porous, 
fractal aggregates. Here, we use the porosity recipe of Ormel 
et al. (2007) , who describe this fractal growth by their en- 
largement parameter as 



*new = 



mp^p + mt^t 



X l-f 



mt^t 



0.425 



+ *add , (47) 



where ^'add is a correction factor in case of rrip « too and 
otherwise zero (for details see their Sect. 2.4). This equation 
predicts an increasing porosity in every hit-and-stick (SI) 
collision. In collisions that lead to sticking through surface 
effects (S2), we assume that the compaction of the aggre- 
gates is so little, that their porosity is unaffected. So the 
aggregates are merged and only the mass and volume of 
both are being added, thus, 



(48) 



One exception for the sticking through surface effects (S2) 
occurs in a small parameter space which is determined by 

the experiments of Blum & Wurm (2000). For the smallest 
masses and a velocity around 100 cm s^^, Blum & Wurm 



(2000) found sticking of fractal aggregates in the 'pP' and 
'pC regimes that goes with a restructuring and, thus, com- 
paction of the projectiles. In this case, we assume a com- 
paction of the projectile by a factor of 1.5 in volume filling 
factor, thus 



Vt^t + min (1.5Vp0p, ^c) 



(49) 



An increasing filling factor is also applied for sticking by 
deep penetration (S3). Here, the mass of the projectile is 
added to the target while the new volume must be less than 
Vt -|- V^. The new volume filhng factor will be 



14< 



pV^p 



(50) 



where Vnew is taken from Eq. 17 (compact projectile) or as 

Kiom' = Vi— Vci. with Vcr. from Eq. 20 (porous projectile). 
In the cases where we transfer mass from one aggregate to 
the other, we always assume that this mass is previously 
compacted by a factor of 1.5 in volume filling factor, but 
cannot be compacted higher than the critical filling factor 
(pc- For the bouncing with mass transfer (B2) we have good 
arguments for this assumption as this compaction is con- 
sistent with XRT measurements of Giittler et al. (2009), 
who also showed that it is likely that this compacted ma- 
terial is transferred to the projectile (see their Figs. 7 and 
9). Without better knowledge, we assume the same com- 
paction of transferred material for fragmentation with mass 
transfer (F3 and S4) and for these three cases we again use 
Eq. 49. Here, we have to note that in the case of bounc- 
ing with mass transfer (B2) and fragmentation with mass 
transfer (F3) the indices of target and projectile need to be 
swapped as the projectile is accreting mass in this coUisional 
outcome. For the fragments in S4 and F3 as well as for those 
in the case of Fl and F2, we assume an unchanged porosity 
with respect to the destroyed aggregate. The most sophis- 
ticated compaction model is used for collisions that lead to 
bouncing with compaction (Bl). Although Weidling et al. 
(2009) measured the compaction only for a small range of 
aggregate sizes and collision velocities, they derived an an- 
alytic model to scale this compaction in collision velocity 
and showed that it is independent in aggregate mass. We 
follow this model but release it from the experimental bias 
due to the (j) = 0.15 samples they used. As outlined in detail 
in Sect. 3, we basically use Eq. 25, and scale the (/'max(^) 
according to Eq. 31 (furthermore using Eqs. 26 30). 

In summary, one can say that the aggregates' porosities 
can only be increased by the coUisional outcomes SI, S4, 
and F3 (see Table 3), where the hit-and-stick (SI) collisions 
will have the most effect. While some coUisional outcomes 
are neutral in terms of porosity evolution (Fl and F2), the 
main processes which lead to more compact aggregates are 
S3 and Bl. 



6. Discussion 

In the previous sections we have developed a comprehensive 
model for the coUisional interaction between protoplanetary 
dust aggregates. The culmination of this effort is Fig. 11, 
which presents a general collision model based on 19 dif- 
ferent dust-collision experiments, which will be adopted in 
Paper II. Since it plays a vital role, it is worth a critical 
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Table 3. Overview on the porosity evolution in the different 
coUisional outcomes. 



collisional outcomes porosity evolution equation 



SI 


fluffier 


47 


S2 


neutral or compaction 


48 or 49 


b3 


compaction 


17, 20, 50 


S4 (target) 


fluffier 


49 


S4 (projectile) 


neutral 




Bl 


compaction 


25 - 31 


B2 (target) 


neutral 




B2 (projectile) 


both 


49*^ 


Fl 


neutral 




F2 


neutral 




F3 (target) 


fluffier 


49*^ 


F3 (projectile) 


neutral 





*'The indices of target and projectile must be swapped here. 



appraisal. In a few examples, we want to discuss the main 
simplifications and shortcomings of our current model. 

(1) The categorization into collisions between similar- 
sized and different-sized dust aggregates (see Figs. 10 and 
11) is well-motivated as we pointed out in Sect. 4. However, 
we may ask ourselves whether this binarization is funda- 
mentally correct, if we need more than two categories, or 
'soft' transitions between the regimes. At this stage, a more 
complex treatment would be impractical due to the lack of 
experiments treating this problem. 

(2) The binary treatment of porosity (i.e. (f> < (pc for 
'porous' and > 0c for 'compact' dust aggregates) is also a 
questionable assumption. Although we see fundamental dif- 
ferences in the collision behavior when we use, e.g., porous 
or compact targets, there might be a smooth transition 
from the more 'porous' to the more 'compact' collisions. 
In addition to that, the assumed value <^c = 0.4 is rea- 
sonable but not empirically affirmed. On top of that, the 
maximum compaction that a dust aggregate can achieve in 
a collision depends on many parameters, such as, e.g., the 
size distribution of the monomer grains (Blum et al. 2006) 
and the ability of the granular material to creep sideways 
inside a dust aggregate (Giittler et al. 2009). 

(3) Although the total number of experiments, upon 
which our model is based, is unsurpassedly large, the total 
coverage of parameter space (see the experiment boxes in 
Fig. 11) is still small. Thus, we sometimes apply extrap- 
olations into extremely remote parameter-space regions. 
Although not quantifiable, it must be c;lear that the er- 
ror of each extrapolation grows with the distance to the 
experimentally confirmed domains (i.e. the boxes in Fig. 
11). Clearly, more experiments are required to fill the pa- 
rameter space, and the identification of the key regions in 
the mass-velocity plane is exactly one of the goals of Paper 
II. 

(4) With such new experiments, performed at the 'hot 
spots' predicted in Paper II. we will not only close gaps in 
our knowledge of the collision physics of dust aggregates 
but will most certainly reveal completely new effects. The 
rather simple 'cc' panel in Fig. 11 as compared to the more 
complex 'pC is due to the fact that there are hardly any 
experiments that back-up the 'cc' regime, whereas in the 
'pC case we have a rather good experimental coverage of 
the parameter space. 

In summary, the sophisticated nature of our collision 
model is both its strength and its weakness. The drawbacks 



of identifying four parameters that shape the collision out- 
come are that rather crude approximations and extrapola- 
tions have to be made. However, to acknowledge the role 
of, e.g., porosity through a binary treatment is still better 
than to not treat this parameter at all. Our new collision 
model represents the first attempt to include all existing 
laboratory experiments (for the material properties of in- 
terest); collisional evolution models can enormously profit 
from this effort. 

6.1. The Bottleneck for Protoplanetary Dust Growth 

In this paper, we have presented the framework and phys- 
ical background for an extended growth simulation. What 

is to be expected from this? Here, we can speculate under 
which conditions growth in PPDs is most favorable. A view 
on Fig. 11 immediately shows that large dust aggregates 
can preferentially grow for realistic collision velocities in the 
'cC and 'pC collision regimes (and to a lesser extent in the 
'pc ' case) , due to fragmentation with mass transfer (S4) . 
For this to happen, a broad mass distribution of protoplan- 
etary dust must be present. This prerequisite for efficient 
growth towards planetesimal sizes has also been suggested 
by Teiser & Wurm (2009b, see their Fig. 11). Agglomeration 
experiments with micrometer-sized dust grains and a stick- 
ing probability of unity (experiments 1 - 3 in Table 1) have 
shown that nature chooses a rather narrow size distribution 
for the initial fractal growth phase. If this changes when the 
physical conditions leave no room for growth under quasi- 
monodisperse conditions, i.e. whether nature is so 'adap- 
tive' and 'target-oriented' to find out that growth can only 
proceed with a wide size distribution, will be the subject of 
Paper II, in which we apply the findings of this paper to a 
collisional evolution model. 

6.2. Influence of the Adopted Material Properties 

The choice of material in our model is 1.5 fjoa diameter silica 

dust as most of the underlying experiments were performed 
with this material. Many experiments (Blum & Wurm 2000; 
Langkowski et al. 2008; Blum & Wurm 2008) showed that 
this material is at least in a qualitative sense representa- 
tive for other silicatic materials - also for irregular grains 
with a broader size distribution. Still, the grain size of the 
dust material may have a quantitative influence on the col- 
lisional outcomes. For example, dust aggregates consisting 
of 0.1 jim are assumed to be stickier and more rigid (Wada 
et al. 2007, 2008, 2009), because the grain size may scale 
the rolling force or breaking energy entering into Eqs. 7 
and 12. However, due to a lack of experiments with smaller 
monomer sizes, we cannot give a scaling for our model for 
smaller monomer sizes at this point. Moreover, organic or 
icy material in the outer regions of PPDs or oxides and sin- 
tered material in the inner regions may have a big impact 
on the collisional outcome, i.e. in enhancing the stickiness 
of the material and thereby potentially opening new growth 
channels. 

As for organic materials, Kouchi et al. (2002) found 
an enhanced sticking of cm-sized bodies covered with a 
1 mm thick layer of organic material at velocities as high 
as 500 cm s~^ and a temperature of ~ 250 K. Also icy ma- 
terials are likely believed to have an enhanced sticking effi- 
ciency compared to silicatic materials. Hatzes et al. (1991) 
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collided 5 cm diameter solid ice spheres, which were cov- 
ered with a 10 - 100 /um thick layer of frost. They found 
sticking for a velocity of 0.03 cm s~^, which is in a regime 
where our model for refractory silicatic material predicts 
bouncing (see 'pp' or 'cc' in Fig. 11). Sintering of porous 
dust aggregate may occur in the inner regions near the cen- 
tral star or - triggered by transient heating events (e.g. 
hghtning, Guttler et al. 2008) - even further out. Ongoing 
studies with sintered dust aggregates (Poppe 2003) show an 
increased material strength (e.g. tensile strength) by an or- 
der of magnitude (C. Giittler & J. Blum, unpublished data). 
This would at least make the material robust against frag- 
mentation processes and qualitatively shift them from the 
porous to the compact regime in our model - without neces- 
sarily being compact. Due to a severe lack on experimental 
data for all these materials, it is necessary and justified to 
restrict our model to silicates at around 1 AU while it is 
to be kept in mind that these examples of rather unknown 
materials might potentially favor growth in other regions 
in PPDs. 
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